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Abstract 

^ New high resolution spectra of the short period (P ~1.76 days) young detached 

binary LT CMa are reported for the first time. By combining the results from the 
analysis of new radial velocity curves and published light curves, we determine values 
for the masses, radii and temperatures as follows: Mi = 5.59 (0.20) M , R\ = 3.56 
(0.07) R Q and T cffl = 17000 (500) K for the primary and M 2 = 3.36 (0.14) M , R 2 = 
2.04 (0.05) R Q and T eff2 = 13140 (800) K for the secondary. Static absorbtion features 
apart from those coming from the close binary components are detected in the several 
j_| ■ spectral regions. If these absorbtion features are from a third star, as the light curve 

solutions support, its radial velocity is measured to be RV3 =70(8) km s _1 . 

The orbit of the binary system is proved to be eccentric (e =0.059) and thus the 
apsidal motion exists. The estimated linear advance in longitude of periastron cor- 
responds to an apsidal motion of U = 69±5 yr for the system. The average internal 
structure constant log£; 2j0 bs=-2.53 of LT CMa is found smaller than its theoretical 
value of log/c2,thco=~2.22 suggesting the stars would have more central concentra- 
tion in mass. The photometric distance of LT CMa (d =535±45 pc) is found to be 
much smaller than the distance of CMa OBI association (1150 pc) which rules out 
membership. A comparison with current stellar evolution models for solar metallicity 
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indicates that LT CMa (35 Myr) is much older than the CMa OBI association (3 
Myr), confirming that LT CMa is not a member of CMa OBI. The kinematical and 
dynamical analysis indicate LT CMa is orbiting the galaxy in a circular orbit and 
belongs to the young thin-disk population. 

Key words: stars: binaries: eclipsing - stars: early-type - stars: evolution - 
stars: fundamental parameters - stars: individual (LT Canis Majoris) 



1. Introduction 

Detailed analysis of early-type double-lined eclipsing binaries (dEBs) in OB associations 
yields not only the precise absolute stellar parameters but also information about the associ- 
ation in which they are embedded. The systemic velocity of dEBs, when combined with their 
distances, allows us to determine the kinematics of the association. The evolutionary stage of 
the dEBs projects the age of the association, and the chemical composition of the components 
of dEBs gives the chemical composition of all member stars of the association under the as- 
sumption that all member stars are formed in the same medium. This valuable information was 
successfully obtained for a number of OB associations by applying modern analysis techniques 
on early- type systems (i.e. among others, rj Mus (Baki§ et al., 2007); DW Car (Southworth 
& Clausen, 2007); AG Per (Gimenez & Clausen, 1994)). In the frame of a project aimed to 
investigate early-type dEBs in the vicinity of OB associations, we observed several dEBs in the 
vicinity of the CMa OBI association. Some are found to be members and others not. 

LT CMa (also HIP 34080, HD 53303) is a relatively bright (V ~7.4 mag), short-period 
(-Porb=l-759535 days), eclipsing binary whose light variations were discovered with the Hipparcos 
satellite (ESA, 1997). To our knowledge, there is no detailed photometric or spectroscopic study 
of the system to date. Claria (1974) reported the Johnson colors of LT CMa as B — V = —0.11 
and U — B = —0.57. Otero (2005) comparing the Hipparcos (ESA, 1997) and ASAS (Pojmanski, 
2002) photometric data interpreted an apsidal motion in the system from the shift of the 
secondary minima. Recently, LT CMa is listed in the catalogue of eclipsing binaries with 
eccentric orbits by Bulut & Demircan (2007). 

In the present study, new high-resolution (R ~ 40000) spectroscopic observations are 
used together with all photometric data available from surveys of LT CMa to reveal the apsidal 
motion rate of the system (section 5) and the absolute stellar parameters of the components 
(section 6.1). A discussion on the membership of LT CMa to the CMa OBI association on 
the basis of evolutionary properties is given in section 6.2. The kinematical properties and the 
population type of LT CMa are determined in section 6.3. Finally, we calculated the internal 
structure constants in section 6.4 for a comparison with theoretical predictions. 



2. Spectroscopic Observations and Data Reduction 

LT CMa has been spectroscopically observed using the Coude-Echelle Spectrograph 
(CES) at the 1.5-m RTT150 telescope of the TUBITAK National Observatory (TUG). 
Presently, CES gives spectra with a resolving power of R ~40000 and is able to provide 68 
spectral orders between 3956-8772 A in a single CCD frame. The detector is a SAO RAS 
lkxlk nitrogen cooled CCD camera. Technical specifications of the CES are described by 
Bikmaev et al. (2005). 

In total, 14 spectra have been collected in five observing runs. The log of the observations 
is given in table 1, where the HJD of the observation corresponds to the mid-exposure times. 
For wavelength calibration, Thorium- Argon lamp spectra were taken in the beginning and at 
the end of each run. Each image was dealt with using IRAF 1 according to normal procedures of 
bias substraction, flat-field division, background-light subtraction and wavelength calibration. 
For the orders where very broad stellar lines such as H a and He I (4471 A) are near the order 
edges, the continuum correction for these orders has been performed by means of dividing the 
order by the normalized continuum function of the previous or next order where no stellar lines 
exist. 

3. Radial Velocities and Spectroscopic Orbit 

The components of the close binary show broad lines. In addition, a careful inspection 
of the spectra at different orbital phases reveals at several locations weak, static absorption 
features that do not participate in the Doppler motions. These static features include H a , H<5 
and Mg II (4481 A) and may be due to a third star that was not detected in the He I lines. Since 
the latter are among the strongest absorption lines in B-type stars, the static features may be 
associated with a cooler third star. The composite spectra of the close binary components and 
the static absorption features at three orbital phases are shown in figure 1. 

The radial velocity (RV) of the static absorption features is 71.2 (7.1) km s" 1 from the 
H Q line and 68.2 (7.7) km s^ 1 from the Mg II (4481 A) line. Errors given in brackets are the 
standard deviations of the measurements which were made by means of fitting Gaussian to the 
central part of the line. Finally, the mean RV of is about 70 (8) km s _1 . 

The strong blending of the broad Balmer series hydrogen lines of the (three) components 
makes them useless for RV measurements. We selected the spectral order 53 for RV study, 
where relatively strong He I (4471 A) and Mg II (4481 A) lines reside. We attempted to measure 
the Doppler shifts of individual spectral lines by fitting a Gaussian to the central part of 
individual lines. However, in some spectra no reliable measurement is possible, either due to 



1 IRAF is distributed by the National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National 
Science Foundation. 



Table 1. Journal of spectroscopic observations. Signal-to-noise (S/N) ratio refers to the continuum near 5700 A. 



ID 


HJD 

(-2450000) 


S/N 


Exposure Time 

(s) 


1 


4786.50846 


50 


3200 


2 


4787.44775 


30 


2500 


3 


4787.51254 


30 


2500 


4 


4787.54492 


50 


3000 


5 


4839.42717 


50 


2400 


6 


4839.45624 


50 


2400 


7 


4839.48497 


50 


2400 


8 


4839.56654 


30 


2400 


9 


4840.29345 


40 


3600 


10 


4840.39241 


40 


3600 


11 


4840.48616 


35 


3600 


12 


4840.52931 


20 


3600 


13 


4841.34322 


55 


2400 


14 


4841.37807 


65 


2400 



the lower signal-to-noise ratio or to the smaller Doppler shifts, or to a combination of both. 
At these phases, radial velocities were not measured and they are left blank in table 2, where 
measured RVs are presented. 

Keplerian spectroscopic orbits have been fitted to the RVs by means of a differential 
corrections technique. The ephemeris for primary and secondary minima were derived by Otero 
(2005) as given in Eqs. 1 and 2. We adopted the orbital period of LT CMa from the primary 
minimum ephemeris (Eq. 1) and fixed it during the orbital solutions. Other parameters such 
as velocity semi-amplitudes of the components (-^1,2), systemic velocity (RVy), longitude of 
periastron (a;) and time of periastron passage (To) were converged. The orbital eccentricity of 
LT CMa is close to zero (see table 4). In view of the low number of available spectra and the 
lack of spectra at quadrature phases, e=0.059 was kept fixed at the value obtained from light 
curve analysis (see table 5), rather than to let it converge to a spurious value. The orbital 
solution is presented in figure 2 and table 3. 

Min I(HJD) = 2448388.537+1.759535 x E, (1) 

Min II(HJD) = 2448389.445 + 1.759503 x E. (2) 

For a check of our spectroscopic orbital solution, the spectral disentangling technique 
was used in the region of He I (4471 A) and Mg II (4481 A) lines. The spectral disentangling 
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Table 2. Radial velocities of LT CMa. Ephemeris given in table 3 was used as zero epoch for phase 

calculation. Values in brackets arc the O — C residuals from the theoretical Kcplcrian orbital fitting. 



Time 


Phase 


RVi 


RV 2 


HJD 


<t> 


(km s- 1 ) 


(km s _1 ) 


2454786.50846 


0.184 


-108.3 (4.4) 


242.5 (2.9) 


2454787.44775 


0.718 


- 


- 


2454787.51254 


0.755 


153.2 (7.1) 


-172.1 (10.0) 


2454787.54492 


0.773 


138.4 (1.1) 


-180.9 (10.7) 


2454839.42717 


0.259 


-109.9 (5.0) 


267.1 (8.4) 


2454839.45624 


0.276 


-116.0 (7.6) 


242.9 (10.6) 


2454839.48497 


0.292 


-115.4 (5.4) 


251.5 (0.6) 


2454839.56654 


0.339 


- 


- 


2454840.29345 


0.752 


- 


- 


2454840.39241 


0.808 


134.8 (7.9) 


-165.1 (6.7) 


2454840.48616 


0.861 


- 


-150.6 (9.9) 


2454840.52931 


0.886 


- 


- 


2454841.34322 


0.348 


- 


- 


2454841.37807 


0.368 


-81.8 (6.4) 


- 



Table 3. The adopted spectroscopic orbital parameters of LT CMa. 



Parameter 



Value 



P (days) 1.759535 (fixed) 

T (HJD-2454000) 786.1848 ± 0.0124 



e 

w sp (rad) 
K x (kms - 1 ) 
K 2 (kms - 1 ) 

<1 

RVy (kms" 1 ) 
mi sin 3 i (M Q ) 
m 2 sin 3 i (M©) 
asini (i?©) 
RMS (kms - 1 ) 



0.059 (fixed) 
3.05 ± 0.36 
132.6 ± 2.9 
220.4 ± 3.1 

0.602 ± 0.020 
23.6 ± 2.0 

4.977 ± 0.070 

2.995 ± 0.065 

12.26 ± 0.20 

7.44 
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Fig. 1. H Q (top) and Mg II (4481 A) (bottom) lines of the components at three orbital phases. In each 
panel, I, II and III represent the lines of the primary, the secondary and the third stars, respectively. 
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Fig. 2. Theoretical curve fitting to the RV data, 
and secondary stars are shown with filled and 
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empty circles, 
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respectively. 



technique has the advantage to find the spectroscopic-orbit elements and the component spectra 
simultaneously without requiring the intermediate step of radial velocity measurements on 
stellar lines (cf. Hadrava, 1995). It yields, in our case, perfectly consistent results with the RVs 
obtained from Gaussian fitting. However, spurious patterns in the reconstructed component 
spectra, especially at the edges, indicate the necessity to combine the spectral orders in an 
accurate way, such that the edges of the reconstructed regions are in the continuum, and, 
possibly, the addition of the third component. The present data sets is also non-optimal in the 
sense that it covers the extreme Doppler shifts, but not their whole range, and that the eclipses 
are not covered. The different line-dilution in mid-eclipse spectra are required to obtain a good 
solution of the low-frequency Fourier components of the component spectra. Summarizing, a 
successful exploitation of the disentanglement techniques will be undertaken when the data set 
is appropriately extended. 

4. Analysis of Light Curves 

4-1. Binary Model and Input Parameters 

Analysis of the light curves was achieved on the basis of three different photometric 
data sets (f/ p -band data of Hipparcos (ESA, 1997), y-band data of ASAS (Pojmanski, 2002) 
and INTEGRAL surveys). Hipparcos and INTEGRAL observations have in total 68 and 480 
photometric measurements, respectively. The usefulness of the Hipparcos data is limited by 
lack of observations in the secondary minimum and the usefulness of the INTEGRAL data 
is limited by insufficient sensitivity. These data sets were used only for the determination of 
periastron longitude. The 315 photometric measurements of ASAS, including eclipses, were 
used for the determination of light-curve elements. All available photometric data are phased 
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Fig. 3. Hipparcos, ASAS and INTEGRAL light curves of LT CMa. 

Observing seasons are indicated with survey name using different symbols. 

and re- scaled in figure 3. 

Otero (2005), using Hipparcos and ASAS photometric data, commented the existence of 
apsidal motion in the system without mentioning its rate. Apsidal motion may lead to wrong 
light curve elements when photometric data is not carefully studied. Especially, when the data 
is spread over years, one should pay attention to analyze subsets separately, for example each 
observing season separately, in order to reduce the influence of apsidal motion in a specific 
set to an insignificant level (cf. Baki§ et al. 2008). In the present paper, we use 5 subsets 
of ASAS data over a time base of 8 years (2001 - 2009). The 2003 and 2004 data covering 
the whole orbital cycle, were found to be more suitable for a reliable light-curve analysis. We 
discarded subsets of 2005, 2006, 2009 due to either low number of data or relatively large 
scatter. Except for the 2007 data sets of ASAS and the INTEGRAL, which were analyzed 
together, we performed a separate analysis for each ASAS and Hipparcos data set using the 
2003 version of the Wilson & Devinney (WD; Wilson & Devinney 1971; Wilson 1994) program. 

The effective temperature of the primary component must be determined for the light- 
curve modeling. The Q-method of Johnson & Morgan (1953) has been used to estimate 
unreddened colours and interstellar reddening. Using U — B = —0.57(0.01) and B — V = 
-0.11(0.009) colours given by Claria (1974), Q = -0.49(0.02) and the unreddened colour 
(B -V) = -0.17(0.03). 
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Consequently the colour excess E(B - V) =0.06(0.03), (U - B) = -0.61(0.03) and the 
visual extinction towards LT CMa Ay =0.19(0.03) mag. These unreddened colours correspond 
to an interpolated spectral type of B4.5V (Fitzgerald, 1970), which corresponds to a temper- 
ature of 16200 K according to the calibration tables of Siess, Forestini & Dougados (1997). 
However, it should be noticed that these colours and the corresponding temperature are ob- 
tained from the combined light of components including the third companion, resulting in a 
slightly redder colour and lower temperature than the intrinsic colour of the primary compo- 
nent. The light curve and the spectral lines show that the light of the primary star dominates 
largely. Therefore, it is reasonable to adopt a temperature of 17000 K which is also consistent 
with the spectral type of the primary (see section 6.1). 

4-2. Determination of the Photometric Elements 

With the analysis strategy of the light curves and the input parameters determined 
in section 4.1, the solution of the light curves has been performed. The temperature of the 
primary was fixed at T e fn=17000 K and the temperature of the secondary was left to converge 
to T e g2- The third light contribution to the total flux of the system was set initially 0.1, based 
on the line strengths of the static component in H a and Mg II (4481 A). Gravity darkening 
exponents g± — gi = 1 and bolometric albedos A\ — Ai — 1 were set for radiative envelopes (von 
Zeipel, 1924). Logarithmic limb-darkening law was used and limb-darkening coefficients were 
taken from van Hamme (1993). The surface potentials (^1,2), light factors of the components 
(Zi j2 ), orbital inclination (i), eccentricity (e) and longitude of periastron (a;) were the adjustable 
parameters during the light curve modeling. The initial guess of the longitude of periastron (u) 
was taken from the spectroscopic solution in section 3. Mass ratio was fixed at the spectroscopic 
value of q = 0.602. 

The solutions of V-band ASAS light curves of 2003 and 2004 years converged very 
rapidly, and had the smallest residuals. The third light contribution did not deviate signifi- 
cantly from the initial guess. Input parameters were then altered to check the consistency and 
uniqueness of the solution. These new input parameters were converged to the parameters of 
the first solution, which shows the consistency of solutions listed in table 4. 

Using the light curve elements from ASAS data sets (2003, 2004) as fixed parameters, 
Hipparcos, INTEGRAL and ASAS (2001, 2007, 2008) light curves were analyzed to determine 
the longitude of the periastron. Since the Hipparcos data are also spread over three years, 
the Hipparcos measurements were also divided into two data sets for the years 1990 and 1992. 
During the fitting to the observational data sets (1990, 1992, 2001, 2003, 2004, 2007 and 2008), 
the longitude of the periastron and the light factors of the components were the parameters 
adjusted. The analysis of each light curve set yielded different values of u (see section 5) which 
clearly showed a time dependent apsidal motion. Apsidal motion is investigated in section 5 in 
more detail. 



Table 4. Light curve fitting results for LT CMa. The parameters converged during the fittings are shown 
with their errors. The parameters without errors are those adopted from the solutions of 2003 and 2004 
light curves (last two columns). The goodness of the fits to each data sets are given in the bottom of 
the table as x 2 together with the number of observations (NOBS) in the respective data set. Note results 
for 2007 are from two data sets ASAS and INTEGRAL. NOBS and \ 2 are given in brackets accordingly. 



Parameters 


Hpi989 


-Hpi991 


V2001 


V2007 


V2008 


V2003 


V2004 


2i(K) 


17000 


17000 


17000 


17000 


17000 


17000 


17000 


T C (K) 


13140 


13140 


13140 


13140 


13140 


13140±210 


13135±315 


q (Af a /Mi) 


0.602 


0.602 


0.602 


0.602 


0.602 


0.602 


0.602 


e 


0.059 


0.059 


0.059 


0.059 


0.059 


0.058±0.003 


0.059±0.006 


w (rad) 


1.23±0.09 


1.39±0.09 


2.40±0.11 


2.80±0.11 


2.60±0.16 


2.39±0.08 


2.55±0.16 


h 


0.757±0.011 


0.744±0.011 


0.737 


0.737 


0.737 


0.732±0.008 


0.743±0.008 


h 


0.163±0.015 


0.161±0.015 


0.162 


0.162 


0.162 


0.164±0.012 


0.159±0.012 


h 


0.080±0.020 


0.095±0.020 


0.102 


0.102 


0.102 


0.104±0.021 


0.098±0.020 


fit 


4.25 


4.25 


4.25 


4.25 


4.25 


4.204±0.041 


4.300±0.055 


(h 


5.06 


5.06 


5.06 


5.06 


5.06 


5.012±0.062 


5.133±0.081 


ri (pole) 


0.274 


0.274 


0.274 


0.274 


0.274 


0.277±0.003 


0.272±0.004 


r\ (point) 


0.289 


0.289 


0.289 


0.289 


0.289 


0.292±0.003 


0.286±0.004 


r\ (side) 


0.280 


0.280 


0.280 


0.280 


0.280 


0.282±0.003 


0.277±0.004 


r\ (back) 


0.286 


0.286 


0.286 


0.286 


0.286 


0.289±0.003 


0.283±0.004 


n (volume) 


0.280 


0.280 


0.280 


0.280 


0.280 


0.283±0.003 


0.277±0.004 


r 2 (pole) 


0.157 


0.157 


0.157 


0.157 


0.157 


0.159±0.002 


0.155±0.003 


r2 (point) 


0.159 


0.159 


0.159 


0.159 


0.159 


0.162±0.002 


0.157±0.003 


ra (side) 


0.157 


0.157 


0.157 


0.157 


0.157 


0.160±0.002 


0.155±0.003 


r-i (back) 


0.159 


0.159 


0.159 


0.159 


0.159 


0.162±0.002 


0.157±0.003 


r-2. (volume) 


0.158 


0.158 


0.158 


0.158 


0.158 


0.160±0.002 


0.156±0.003 


*(°) 


74.2 


74.2 


74.2 


74.2 


74.2 


73.9±0.3 


74.5±0.4 


NOBS 


32 


36 


39 


480, 64 


78 


66 


68 


x 2 


0.001527 


0.002589 


0.003311 


0.147187, 0.011011 


0.011290 


0.001600 


0.002633 



Adopted light curve solutions for each data sets of LT CMa are shown in figure 4. The 
light curve elements obtained from these data sets are given in table 4. For clarity, we give the 
adopted light curve elements in table 5. 

5. Apsidal Motion 

The orbit of LT CMa has been determined to be slightly eccentric (e = 0.059) from the 
solutions of the light curves in section 4. This indicates that LT CMa is a system showing 
apsidal motion. However, due to the small number and unfavorable distribution of the times 
of minima, it is not possible to determine the apsidal motion period by the O — C variations. 
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Table 5. Adopted V-band light curve elements for LT CMa. f! C r stands for the potential of the inner critical Roche lobe. 



Parameter 


Primary Secondary 


T(K) 


17000 13140±315 


q(M 2 /M 1 ) 


0.602 


i(°) 


74.2±0.3 


e 


0.059±0.006 



/(/i, 2 //Totai) 0.737±0.008 0.162±0.012 

/(/ 3 //Total) 0.102±0.020 

fl 4.25±0.04 5.06±0.06 

cr •. ) . _ _ 

r 0.280±0.004 0.160±0.003 



n i ' i ' i ' r 



tfXi^ \ f> — ■ U " "" " J ' a" " " 
ASAS - 2008 



"i ' r 




•%>r 



#> — • * — *-< 



'o 5; <r-*gts «#. 




u e 6^ e r W l u a '^ i 



a an n 9 U n m — % n " a u°u uuu n «» ' u a 



ua" u « 'u^ n 



— »-a — °a 5—0 — 8 5 tjl u 



I n^jn u »n [g— 




0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 

Phase 

Fig. 4. Photometric solutions of V and iJ p -band light curves. 
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Fig. 5. Apsidal motion calculation through a linear regression to longitude of periastron 
(a;) values obtained from the analysis of eight light curves and one RV curve. Least 

squares fit to the u values with their uncertainties yielded the slope of the straight 
line as 0.00025±0. 00002 rad d" 1 . The ±l-u levels is also shown with dashed lines. 

Therefore, determination of the apsidal motion rate has been performed by means of linear 
least-squares fitting to the longitude of periastron (u) values which were obtained from nine 
different data sets (two Hipparcos light curves, five ASAS light curves, one Integral light curve 
and the spectroscopic data). The best linear fit (u = o>o + Co T) to the set of longitude of 
periastron values is presented in figure 5. The data are consistent with a linear advance in 
longitude of periastron, d>=0.0252 ± 0.0020 deg cycle" 1 . This apsidal motion rate corresponds 
to an apsidal motion period of U = 69 ± 5 yr which is a relatively short apsidal period compared 
to the eccentric binary systems given by (Bulut & Demircan 2007). 

6. Discussion 

6.1. Absolute dimensions and distance of LT CMa 

The fundamental astrophysical parameters of LT CMa, which were derived from spec- 
troscopic results in table 3 and from the light-curve results in table 4, are summarized in table 
6. The temperature T e si = 17000 K, mass Mi = 5.59M and radii R\ = 3.56f? of the primary 
component correspond to a spectral type B4V. The temperature T eff2 = 13140 K of the sec- 
ondary star implies a B6-B7 spectral type ZAMS star which is also very consistent with its 
mass of M 2 = 3.36M and radius of R2 = 2.O4i? (i.e. Straizys & Kuriliene 1981). 
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The unreddened Johnson ^-magnitude (Clarie 1974) of LT CMa, when combined with 
the light contributions as derived from the light-curve analysis, including the distant companion, 
yields the intrinsic ^-magnitudes myi =7 m .55, mv2 =9 m .19 and mv3 =9 m .72, respectively. 
Using M v =4 m .75 as the absolute visual magnitude of the Sun and bolometric corrections 
BC\ =-1.58 and BC2 =-0.93 mag for the primary and the secondary, from Straizys & Kuriliene 
(1981), bolometric and absolute visual magnitudes of the close binary components are derived 
(see table 6). The visual magnitude and distance modulus indicate a photometric distance of 
535 ± 45 pc to the LT CMa, which is smaller than the distance (826±326 pc) derived for LT 
CMa by van Leeuwen (2007) from the re-analysis of raw Hipparcos data, but not inconsistent 
with it in view of the large error bar on the Hipparcos distance. In the original Hipparcos 
catalogue, the distance to LT CMa is given to be 591±386 pc. Among these measurements, 
the photometrically determined distance of LT CMa comes out to be more precise and reliable. 

The distance of the CMa OBI association (1150 pc) (Claria 1974) is much greater than 
the distance of LT CMa (535 ± 45 pc), which weakens the probability that LT CMa belongs 
to the association. 

The projected rotational velocities of the components were derived to be Hoti x sin(i) 
= 105(10) km s _1 and V^. ot2 xsin(i) = 65(10) km s" 1 by fitting synthetic spectra using Kurucz 
(1993) atmosphere models and the ATLAS9 code (Kurucz 1993). Using the orbital inclination 
(z=74.2 deg) of LT CMa, the observed rotational velocities of the components are found to be 
Vroti = 109(10) km s _1 and V TOt 2 = 67(10) km s _1 which are in agreement with the synchronous 
rotational velocities within the uncertainties. 

Although the light curve solution supports an extra light in the system, the assumption 
of the third star is still preliminary due to the low number of high S/N spectra used in present 
study. Thus, it is desired in the future study to obtain high S/N spectra evenly distributed over 
the whole orbital cycle in order to reveal the identity of additional features in the composite 
spectrum. 

6.2. Evolutionary stage and age of LT CMa 

We have investigated the evolutionary status of LT CMa in the plane of log T e g - logg 
and Mass - logT c g (figure 5) using the latest evolutionary models and isochrones of Girardi et al. 
(2000), which include mass loss and moderate overshooting. Assuming a solar metal content, 
we prepared a set of isochrones corresponding to Y = 0.28 and Z = 0.02. From figure 5a, it can 
be interpreted that both component are consistent with the evolutionary tracks calculated for 
their masses. The primary star is evolved from the ZAMS, whereas the secondary is still on the 
ZAMS. The isochrones of 30 Myr and 40 Myr shown in figure 5a and figure 5b implied a mean 
age of 35(5) Myr for the system. The age of LT CMa derived here is much greater than that 
derived by Claria (1974) for CMa OBI (3 Myr), which rules out membership, when combined 
with the distance difference of LT CMa and the association. 
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Table 6. Close binary stellar parameters of LT CMa. Errors of parameters are given in parenthesis. 



Parameter 


Symbol 


Primary 


Secondary 


Spectral type 


Sp 


B4 V 


B6.5 V 


Mass (M Q ) 


M 


5.59(0.20) 


3.36(0.14) 


Radius (Rq) 


i? 


3.56(0.07) 


2.04(0.05) 


Separation (Rq) 


a 




12.7(0.2) 


Orbital period (days) 


P 




1.759535 


Orbital inclination (°) 


i 




74.2(0.3) 


Mass ratio 


'I 




0.602(0.020) 


Eccentricity 


e 




0.059(0.006) 


Surface gravity (cgs) 


logs 


4.081(0.014) 


4.347(0.019) 


Integrated visual magnitude (mag) 


V 


7.42(0.01) 




Integrated colour index (mag) 


B-V 


-0.11(0.009) 




Colour excess (mag) 


E(B~V) 




0.06(0.03) 


Visual absorption (mag) 


A v 




0.19(0.03) 


Intrinsic colour index (mag) 


[B-V) 




-0.17(0.03) 


Temperature (K) 


T c ff 


17000(500) 


13140(800) 


Luminosity (Lq) 


log L 


2.98(0.03) 


2.05(0.07) 


Bolometric magnitude (mag) 


M bo i 


-2.70(0.05) 


-0.36(0.08) 


Absolute visual magnitude (mag) 


M, 


-1.12(0.05) 


0.57(0.07) 


Bolometric correction (mag) 


BC 


-1.58 


-0.93 


Velocity amplitudes (km s _1 ) 


Kl,2 


132.6(2.9) 


220.4(3.1) 


Systemic velocity (km s~ ) 


RV 1 




23.6(2.0) 


Computed synchronization velocities (km s^ 1 ) 


V synch 


102 


59 


Observed rotational velocities (km s _1 ) 


Kot 


109(10) 


67(10) 


Distance (pc) 


d 




535(45) 


Proper motion (mas yr _1 ) 


flaCOsS, [IS 


-6.82(0.83), 3.53(0.71)* 


Space velocities (km s _1 ) 


u,v,w 


-27.46(2.12), 


-5.21(2.15), -12.45(2.24) 



from Hipparcos catalogue (van Leeuwen 2007). 
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Fig. 6. Locations of the components stars in the plane of logl^g- - logg (upper panel) and in 
the Mass - logT c g (lower panel). Primary and secondary stars are shown with filled and empty 
circles, respectively, together with error bars. Isochrones are for 10 Myr (dotted line), 30 Myr 
(dash line) and 40 Myr (dotted dash line). Evolutionary tracks for the primary and the sec- 
ondary stars are presented with solid lines where the dashed lines around evolutionary tracks are 
due to the uncertainties in masses of the components. The horizontal solid line is for ZAMS. 
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6.3. Kinematical analysis and population type of LT CMa 

To study the kinematical properties of LT CMa, we used the system's centre-of-mass 
velocity, distance and proper motion values, which are given in table 6. The proper motion 
data were taken from the newly reduced Hipparcos catalogue (van Leeuwen 2007), whereas the 
center-of-mass velocity and distance are obtained in the present study. The system's space ve- 
locity was calculated using the Johnson & Soderblom's (1987) algorithm. The u, v and w space 
velocity components and their errors are listed in table 6. To obtain the space velocity precisely, 
the first-order galactic differential rotation correction (DRC) was taken into account (Mihalas 
& Binney, 1981), and -10.33 and -0.90 km s _1 DRCs were applied to u and v space velocity 
components, respectively. The w velocity is not affected in this first-order approximation. As 
for the local standard of rest correction, Mihalas & Binney's (1981) values (9, 12, 7) km s _1 
were used and the total space velocity of LT CMa was obtained as S = 12.45(3.76) km s" 1 . To 
determine the population type of LT CMa, the galactic orbit of the system was examined. Using 
the N-body code of Dinescu, Girardi & van Altena (1999), the system's apogalactic (R ma x) and 
perigalactic (R m in) distances were obtained as 9.61 and 8.35 kpc, respectively. Also, the max- 
imum possible vertical separation of the system from the galactic plane is \z max \ = \z min \ = 30 
pc. The following formulae were used to derive the planar and vertical ellipticities: 

(4) 

The planar and vertical ellipticities were calculated as e p = 0.07 and e v = 0. These values 
show that LT CMa is orbiting the Galaxy in a circular orbit and the system belongs to the 
young thin-disc population. 

6.4- Internal Structure 

The internal structure constant, (fe), is a measure of the density concentration of the 
components of eclipsing binary stars with eccentric orbits. It is an important parameter of 
stellar evolution models. If the orbital eccentricity, the orbital and apsidal motion periods, the 
masses and radii of the components are known, the mean internal structure constant of the 
component stars in the system can be derived. 

The observational mean value of the internal structure constant, (&2,obs), is given by the 
following formula (Kopal 1978): 

k2,obs = ; 77, (5) 

c 2 i + c 22 U 

where P and U are the orbital and apsidal period, and the constants c^ (i= 1, 2): 

<* = r t [(— ) (1 + — )/(e) + -jf-9(e)}, (6) 
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6p 


-K'max ~r ±i"min 


e v — 


\ Z-max \ ' ~ \ ^min \ 


-K-max ~T J^-min 



where 

fie) = - 

1 -e 



/(e) = 7^5yP (7) 



(8 + 12e 2 + e 4 )/(e) 2 - 5 
s( e ) = g • (8) 

Here r ; , M- x and u; r .i are the relative radii, masses and the angular velocities of the axial rotation 
of the components, respectively, and u^ is the mean angular velocity of orbital motion. The 
first term in Eq. 6 represents the contribution of rotational distortions, and the second term 
corresponds to the tidal effects. 

It is well known that the observed /c2,obs contains not only the Newtonian contribution 
to apsidal motion, which is due to the stars being distorted extended shapes, but also a small 
amount of relativistic apsidal motion. The theory of General Relativity estimates the relativistic 
contribution to the observed rate as the following Einstein formula (Gimenez, 1985): 



1 M 1 + M 2 , 2/3 
1 - e 2 1 P 



u rel = 5.45 x 10" 4 - -A iZ 2 ) 2 / 3 (9) 



where Mj denotes the masses of the components in solar mass and P is the period of the orbit 
in days. 

The individual internal structure constants, k 2 i, must be combined using the equation: 

j C2lk 2 l t theo + C22k22,theo , 1rA 

K2,theo — j • (,-WJ 

C21 + C22 

to find the weighted average coefficient which is directly comparable to the observed value. 

The observed apsidal motion period of LT CMa U = 69 ± 5 yr, corresponding to a 
total rate of u = 0.0252 ± 0.0020 deg cycle -1 was obtained in section 5. The relativistic 
contribution to the apsidal motion calculated from general relativity is 0.00162 deg cycle -1 or 
amounting to about 6% of the observed rate. After correcting for this effect, an average internal 
structure constant of LT CMa is derived, with assumption of the periastron-synchronization, 
to be logk2,obs=— 2.53. 

The theoretical internal structure constants for the components (log &2i,theo) are taken 
from the new evolutionary models of Claret (2004) with the standard chemical composition 
of (X, Z) = (0.70, 0.02). Among the tabulated values, interpolation for the masses and logg 
of the component stars of LT CMa yielded -2.23 and -2.18 for the primary and the secondary 
component, respectively. The average theoretical value of internal structure constant was 
computed as /c2,theo=-2.22. This value seems significantly different from the observed value, in 
the sense that the component stars appear to be more concentrated in mass than theoretically 
predicted from the evolutionary models. There are several possible reasons for disagreement 
in the case of LT CMa. Perhaps the most important is the fact that the range of obser- 
vations used to derive the apsidal motion period is small compared to the apsidal motion period. 
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